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Abstract 

We obtain a variety of series and integral representations of the digamma 
function tp{a). These in turn provide representations of the evaluations ip{p/q) 
at rational argument and for the polygamma function The approach is 

through a limit definition of the zeroth Stieltjes constant 70(a) = —■ip{a). Sev- 
eral other results are obtained, including product representations for exp[7o(a)] 
and for the Gamma function T{a). In addition, we present series represen- 
tations in terms of trigonometric integrals Ci and Si for ^l^{a) and the Euler 
constant 7 = —ip{l). 
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Introduction and statement of results 

We recall the defining Laurent expansion of the Stieltjes constants 7fc(a)[6l [71 El 

El EDI 122] 

C(^, «) = — T + ^ ^7n(a)(s - 1)", (1.1) 

* n=0 

where C{^^(^) is the Hurwitz zeta function. Notationally, we let ({s) = Ci^A) be 
the Riemann zeta function [TTl [TSl HZJ |2I], T the Gamma function, ip = F'/F be 
the digamma function (e.g., [T]) with 7 = — '?/'(l) the Euler constant, '?/'*^^^ be the 
polygamma functions [1] , and pFq be the generalized hypergeometric function [2] . 
The Stieltjes constants may be expressed through the limit relation [3] 



(-1)" 
7n(a) = j— hm 



^ ln"(fc + a) ln"+^(Ar + a) 



E 



k + a n + 1 



n>0. (1.2) 



Here, a ^ {0, —1, —2, . . .}. For an asymptotic expression for these constants, even 
valid for moderate values of n, |16j (Section 2) may be consulted. 

In this paper, we obtain various representations of the digamma function via the 
connection 70(a) = —ip{a) [22j. These in turn lead to many special cases, including 
the values ip{p/q) for rational argument, and further imply representations of the 
polygamma functions. We obtain product representations for exp[7o(a)] and F(a). 
We present series representations for InF(a), ipia), and 7 using the trigonometric 
integrals Si and Ci. In addition, we provide several summations over parameterized 
values of Ci. 

The following is an example of representations that we develop. 
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Proposition 1. Let Re a > 0. Then we have 



7o(a) = -V-la) = - / -— — — — -3F2(l,2,2;3,t + a + 2;l)rft-lna (1.3) 
2 Jo [t + a){t + a + 1) 

OO 1 k /1\ 

= -lna-E^E(-ir ' + (1-4) 

= - In a - C u"'^ i— du (1.5) 

JO \u — 1 muj 

= -\na- f C If ^^"''11 7 "^^/ xdi/ (1.6a) 
JO Jo (1 — xy) ui[xy) 

— In a — / tt:-? — ; -^dt. (1.6c) 

/„ f^2 lt+k+a~l\ ^ ' 

\ k ) 







k=2 

In (1.6b), T is the triangle with vertices at (1, 0), (0, 1), and (1, 1). 

Thus from (1.4) we have 
Corollary 1. We have for Re a > 

OO k 

a 



1 CO fc r 



fe+i 



(1.7) 



k=l t=0 

Proposition 1 and this Corollary subsume expressions for 7 given in [18j . 

In addition, from Proposition 1 follow representations for the polygamma func- 
tions, and these include 
Corollary 2. We have for Re a > 



k=i ' i=0 



+ y'^-l ( _1 J_\\^3udu (1.9) 

^ ' ai Jo Vm-1 InuJ ^ ' 
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= {-ly-'^-^^ + /' /' ^^^^P^ln^-^M dxdy (1.10a) 
Jo Jo (1 — xy) 

. (-ir ■ <i-i^ + - .n^-(l - Y),X,Y. (1.10.) 

Hence we obtain representations at positive integer arguments for harmonic num- 
bers Hn = Ylk=i generalized harmonic numbers H^^ = J2^=i l/^*^- For we 
have Hn = ip{n + 1) — = ipin + 1) + 7 and 

Hi^^ = [v.('-^)(n + 1) -^(-i)(l)] , (1.11) 

where ip^^'-^^l) = (-l)^(r - l)!C(r). 

Further, if we integrate — /"7o(6)(i6 = Inr(a) we have 
Corollary 3. We have for Re a > 

00 1 k /T,\ 

In r(a) = a(ln a - 1) + 1 + E ^ E(-l)' U [(^ + W + a) - {i + 1) ln(£ + 1)] 

(1.12) 

= a(lna - 1) + 1 + [\u''~' - 1) ( ^— - p- (1.13a) 

JO \n — 1 mu) mu 

= a na-1 -- na + - n 27r + / « + - ; 1.136 

^ '2 2 ^ ' }o Vm-1 Inw 2/lnM ^ ' 

= a(lna - 1) + 1 + [(^^)" ~ '^^^]~ ''^ dxdy (1.14a) 

JO JO (1 — xy) In (a;?/) 

= a(lna - 1) + 1 + / / ~ ~ ^^ c^Xdr. (1.146) 

^ ^ i xrin^(i-r) ^ ^ 

By inspection, we see that the right sides of (1.12)-(1.14) properly vanish at a = 1. 
The term a(lna — 1) on the right sides there is the leading asymptotic form of Inr(a) 
when a — )■ 00. In writing (1.13b), we have used (e.g., [T9], p. 16) 

1 1 1\ du 1 , . X W"-^ - 1 , 

1 ln(27r), / — au = lna. Re a > 0, 



\u — \ In-u 2/ In-u 2 Jo In-u 

(1.15) 



thereby recovering Binet's first expression for Inr(a). 

Tlien we liave from (1-12) 
Corollary 4. We have for Re a > 



oo k 

r(a) = e'-'^a'^ J] H 

k=ll=0 



e(k\ 



fc+i 



:i.i6) 



{i + iy+\ 

A multitude of summation formulas for the digamma function is known (e.g., 
Section 55, [lOj). Our integral representations in particular well provide a basis for 
further developing such summations. We simply give an example here. In light of the 
asymptotic relation as a — )• oo, 



we may consider the following. 

Corollary 5. For Re a > and Re (« + /?)> 0, we have (a) 



:i.i7) 



E 

n=l 



ip{an + (3) — ln(an + (3) + 



1 



+ 



1 



1 yP/a 
1 — V 



2{an + (3) 12(an + 
1 1 1 



a 



int, 2a 12a' 



liav 



dv, 



:i.i8) 



and (b) for a ^ 1 



E 



n=l 



ip{an + (3) — ln(ar;, + (3) + 



1 



k=l 



a 



2fc-l/2 2k 



a , 



2{an + I3)_ 

where Bj is the jth Bernoulli number. 

Then we may determine the asymptotic dependence to all orders of a certain 
second moment of the Riemann xi function. For this, we put ^(s) = (s — l)7r~'^/^r(l + 



s/2)C(s), and E{t) = ^(1/2 + it). We have 
Corollary 6. We have as a — oo 







^ ft\^ fit-l\ 2 cos[(t/2) In a] 



r 



4 



t2 + 1 



dt 



^3/2 oo 1 n 

._|,„„ + ,„p,)_-,, + .»/'^|;__|c(2.,. 

The moment integral here, going back to Ramanujan, is of interest from many points 
of view [1]. 
Let 

ci(^) =- r ^dt, si(^) = r ^-^dt. (1.19) 

Jz t Jo t 

Then we have 
Proposition 2. We have 

1 lf°°l 001 ^ 

7 = - + Ci(vr) + - E — tI^ - 2Si(vrj)] + E - 2Si(vrj)] • (1.20) 

In this expression, Ci(7r) ^ 0.07366079 and therefore the sum terms provide small 
corrections. 

Proposition 3. We have for Re a > 

oo 

2p{a) = Ina + E [2 cos(27rja)Ci(27rja) - sin(27rja)[7r - 2Si(27r ja)]] . (1.21) 

2a .^^ 



Then from %p{n + 1) = ip{n) + 1/n we have 
Corollary 7. We have 



2 oo 

Hn = lnn + 7 + — + 2ECi(27rnj). 
2n 



Alternatively, this result follows from the Euler-Maclaurin summation expression 



i/n = lnn + 7 + — + / — V^rfx = lnn + 7 + — + - / 



1 1 f°° PAnv) , 
-^—dv, 



where Pi is the polynomial given in (2.29). Similarly, for r > 1 we have 



1 n^~'' 

= C(r) + — - + r 

2n^ r — 1 



^1 ^ tp poo 



X 



r+1 



+ — 



Pi{nv) 



dv. 



Let Bn{x) be the Bernoulli polynomial of degree n, Bn = -Bn(O) the nth Bernoulli 
number, and -B„(l/2) = (2^"" — l)Bn [Ij (p. 805). Then we have the following. 
Proposition 4. Let < (3 < 27i. Then we have (a) 



n=l 



,Ci(/3n 



n=l 



^(ln2-7-ln/3)C(2) + ^C'(2) + ^. 



(b) 



^ Ci{f3n) 



n=l 



(-,+,„^)c(4)-c'(4)+g(4+^^-^; 



n=l 

(c) for integers k > 1 



192 8 



£ ^ = (7 + ln/3)C(2fc) - C'(2A:) + ^-J^i2n^'' ^ 



71=1 



n 



2{2k)\ 



f 



(1.22a) 
(1.226) 

(1.22c) 



£(-1)"^^ = (ln2 - 77 - 71n/3)^ + ^C'(4) - ^ + (1.22rf) 



B2k ~ ^2fc , 

(1.23a) 

E(-1)"^S^ = (7 + ln/3)(2i-''= - 1)C(2A;) + (1 - 2'-'')C{2k) + 2'-'\\n2)a2k) 



n=l 



n 



2(2A;)! 



-(2vr) 



f 



5^ 



2fc 



' (3v + vr' 
27r 



5^ 



2k 



dv. 



(1.236) 



(d) for Re a > 1, 



~ Ci(/3n) , , , . \ (27r 



- sec — 

a) V 2 ; 



2^ 



C 1 - a, 



(1.24a) 



2^ 



E(-l)"^^ = (7 + ln/3)(2^-" - l)C(a) + (1 - 2i-'^)C'(a) + 2^-\{a) ln2 



+ 



n=l 
1 



(27r)'' /vra 
sec ' 



C 1 - + C 1 - «, 



27r 



27r 



[4r(a) 
e) for Re a > 1, 

OO -1 

E .^^^..„ Ci(/jn) = (7 + ln/3)[(l - 2-)C(a) - 1] + 5(a 



+ (1 - 2^-'^)C(a) 



V 



(1.246) 



+ 



l4r(a) 



csc(7ra) 



— sm 



. ( fiv — T:a 
+ sm 



C 1 - a, 



2 



(2- - l)C(a) -, 



dv 

V 

(1.25a) 



+ 



Jo 



+ cos 



where, 



S(2ri + 1)" 

I4r(a) 

/Sv — 7ra 



CSC 7ra 



Ci(/5n) = (7 + ln/3)[4-'^(C(a, 1/4) - ({a, 3/4)) - 1] + T(a 

/3f + 7ra 



cos 



47r y 



An J 



+ 4-»[C(a,3/4)-C(a,l/4)]|^, 



^ ' t2 (2n + 1)^ 



n=2 



(-l)"lnn 
(2n + 1)« ' 



(1.256) 



(1.26) 
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and (f) for \z\ < 1, 



y 

^in(n + l) 



Ci(^n) = (7 + ln^) 



1 - ln(l - ^) + - ln(l - ;2 



^^in(n+l) 



H — / — 1) ln(l — z)] — -iz — cos /Jv) ln(l — 22;cos j3v + — sin j3v tan~^ 

z Jo 2 Vl"-^ 



zcos^ 



(1.27) 

Proof of Propositions 

We let ■0' be the trigamma function and = V{h-\-n)/T{h) be the Pochhammer 
symbol. 

Proposition 1. Preliminary relations are contained in 
Lemma 1. We have 



oo . 

70 (a) = E /, 



1 



1 



poo 



{t + k + ay {t + k + a + l){t + k + a) 
1 dt 



dt — In a 



(t + k + ay (t + k + a + 1) 



In a 



k=0 



Jo L 



t + a 



dt — In a — —ip{a). 



(2.1) 



Proof. From relation (1.2) we have 



N 



7o(«) = Jim^E 



fe=0 



1 



N 

lim 

^7■ — ^ 



N^oo 



k=0 



k + a 

k + a + 1 



In (TV + a) 



k + a 



In 



oo „c 



k + a 
1 



In a 



(t + A; + a)2 {t + k + a + l){t + k + a) 



dt — In a. 



(2.2) 
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The rest of the Lemma follows easily upon noting (e.g., p. 259 of [Ij) 

OO 



OO -| 

^'(z) = C(2, z) = E T--^^- (2-3) 



We now write from the Lemma 

7o(a) = / —- -> (777 ; -dt-\na. (2.4) 

'^^ ^ Jo {t + a)^t + a+l) ^^{t + k + a)^{t + k + a + l) ^ ^ 

The integrand of (2.2) being absolutely convergent, the interchange of summation and 
integration is justified. In order to achieve hypergeometric form, we note the ratios 

it + a)k t + a {t + a)k {t + a){t + a + l) 



{t + a + l)k t + a + k' (t + a + 2)fc {t + a + k){t + a + k + 1)' 

Upon using the series definition of the function 3F2, we therefore obtain 



(2.5) 



roD I I 

7o(a) = / TT, — ^TT, — T 3^2(1, t + a,t + a;t + a + l,t + a + 2;l)dt-\iaa 

Jo [t + ay (r + a + Ij 

= 9/ 77— ^77— — TT3i^2(l,2,2;3,t + a + 2;l)f/t-lna. (2.6) 
2 JO (t + ay (t + a + 1) 

Here, we have used the transformation [18j (6), valid for Re s > and Re f — t > 0, 

3^2(1, s,t-s + 1, V- 1) = 3F2{l,v-t,v-s-v-t + 1, V- 1). (2.7) 

V — t 

We have obtained (1.3). 
Next we have 

3^2(1, 2, 2; 3, t + a + 2; 1) = ;E ,,V^ , = 2 ^ + 



j^o (3).(i + a + 2), f^^ U + 2){t + a + 2)j 

of- 1 j! _oy^ 1 j!(t + a + j + l) 

^,j + l(t + a + 2),_i hj + l it + a + 2), " ^"^ 
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The integral in (2.6) becomes 

1 3i^2(l,2,2;3,t + a + 2;l) dt _ f°° ^ j\ {t + a + j + l)dt 



2 Jo (t + ay (t + a + l) Jo {j + 1) {t + a){t + a + l){t + a + 2)j 

(2.9) 

Therefore, from (2.6) we find 

^(a) = -\na+ Y- ■ ^ ; -dt. (2.10) 

Jo ^^{k + l){t + a){t + a+l)---{t + k + a) ^ ' 

We can carry out the integration by using the partial fraction decomposition [13] 



N\ ^ (N\ (-1)^ 

-EL (2-11) 



x[x + 1) ■ ■ -{x + N) f^Q\k){x + k)' 
We then have (1.4). 

If we employ the Beta function integral 

/ u'^'^-\l - ufdu = B{t + a,k + l)= ■ , (2.12) 

^0 [t + aj{t + a + Ijk 

in (2.10) we find 

7o(a) = -lna+/ Y '-dudt (2.13) 

"'0 Jo k + 1 

= -\na + r r M^+^'-i ( - 1 1 dudt. (2.14) 
Jo Jo \u — 1 J 

Performing the integral over t gives (1.5). 
We may rewrite (2.13) as 



1 rl-u y^t+a-l y 



7o('^) = — lna+ / / / 7 TT -dvdudt 

Jo Jo Jo (1 — n) (1 — V) 

= -lna+ / . (2.15a) 

Jo Jo [u — 1) [1 — V) Inu 

11 



We now put u — xy, v — 1 — x and use the Jacobian d{u,v)/d{x,y) — x to obtain 
(1.6a). 

Putting^ = X andy = {1—Y)/X, with Jacobian of transformation d{X,Y) /d{x,y) 
-X, in (1.6a) yields (1.6b). Using -l/ln(l -Y) = - Yfdt in (1.6b) gives 

7o(a) = -^/;(a) = -lna+ / / / ^ {I - X)dtdXdY 



It Jo XY 

roc rl rl (1 - X) 



-\na + r [ [ ^ ^(1 - Yf+^'-^dXdYdt 

Jo Jo Ji-Y XY 

roo fx _ y^+ft-l 



poo rl I Y '-"r"— 

= -lna-/ / ^ [Y + ln(l-Y)]dYdt (2.15b) 

Jo Jo Y 

POO °° 1 rl 

= -lna-/ E7:/(l- Yf+'^-^Y'^-'dYdt 

'fO ^ 2 '^Q 

poo ^ I 

= - In a - / y yBik, t + a)dt, 
where we employed the Beta function integral of (2.12). Then we have the expressions 







roo - Y{k) dt 
= -Ina- / 2^ 



.2 ^ it + a)k 



that are equivalent to (1.6c). 

Corollary 5. From (1.5) we have 



^(a) = Ina - - ttt^ + C u"'^ ( ^— - + ^ - ^Iriu) du. (2.16) 
^ 2a 12a2 io \u-l \nu 2 12 J ^ ^ 



We then interchange summation and integration to find 



E 

n=l 



t/'(an + /3) - ln(an + /3) + — — — — + ^ 



2{an + /3) 12(Q;n + /3)2 
12 



1 ^a+p-l 
1 - 



1 



1 



1 



1 



du. 



(2.17) 



u-l liiM 2 12 
Making the change of variable f = gives part (a) of the Corollary. For part (b), 
we use the asymptotic relation for the digamma function as z — > oo 



In z 



E 



2k 



2z f^2fcz2fc' 
fe=i 



Corollary 6. This follows from the /3 = case of Corollary 5(b), together with the 
reciprocity relation with a, /3 > and a/? = 1 [1] 



7 — ln(27ra) 



2a 



71=1 



7 - ln(27r/3) 



1 



2/3 

2 cos[(t/2) Ina^ 



t2 + l 



-(it, 



TT^/^ JO 

where 0(x) = ?/'(x) + l/2a; — Inx. 

An alternative approach is to expand the integrand of Corollary 5(a), using 



,1/q 



followed by the use of the integrals /o^pn-' v/ (1 — v)]dv = {—lyjKU + !)• 

We note that the reciprocity relation itself provides the complementary asymptotic 

relation as a — 0. For then /3 = 1/a — )■ oo. 

We also note that the leading term of the asymptotic relation in Corollary 6 is 

connected with the skew self-reciprocal inverse Fourier cosine transform 



-1 



In a' 



TT 



._L(,„, + -, + _ + 21n2 
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This transform may be calculated by logarithmic differentiation with respect to x of 
the integral Jl^ cos (at) da with — 1 < Re x < 0, at x = —1/2. 
Remarks. If in (1.5) we put u = exp(— t), we have 



oo / 1 1 

at I ^ 



7o(a) = -Ina + y^ e""* (l^^ " ij ^^'^^^ 

thereby recovering formula 8.361.8 of [T2] (p. 943) for the digamma function. This 
formula is also recovered if we use an integral representation for In, In^; = /o°°(e~* — 
e~^^){dt/t) for Re z > 0, in (1.4). For then we have 



If 



°o 1 poo p-a* 



t 



a - e'Tdt 



= - In a - /°° e-"* (- ^—-] dt. (2. 19) 

JO \t 1 — e~v 

If in (2.15b) we instead carry out the integration over Y, we recover the integral 
of (2.1). 

Special values of 70 are 7o(l/2) = 7 + 2 ln2, and 



1\ ""-^ 1 

+ -)=, + 21„2-2i;— . (2.20) 

Further, 7o(l/4) = 7 + 7r/2 + 3 In 2, and 7(3/4) = 7 — 7r/2 + 3 In 2, giving for instance 
7o(l/4) —70(3/4) = TT. More generally, for rational arguments we have the following 

m, p. 13 or m, P 19). 
Theorem. 

= -7 - - cot — - In g + 2 y cos In f 2 sin —\ , (2.21) 

2 g q \ q J 
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where < p < g; X)' means that when q is even the term with index n = q/2 is 
divided by 2. 

Therefore, we have found representations for all the values of (2.21). 

The digamma and polygamma functions satisfy many properties including func- 
tional equations, duplication and multiplication formulas, and reflection formulas. All 
such properties must be inherent in our various series and integral representations. 
As an illustration, we have 

Corollary 8 (multiplication formula). For integers m > 1 we have 



7o (ma) 

Proof. From (1.5) we have 



1 '"^^^ ( k\ 
lnm + -V7o a + - . (2.22) 

m ^0 V 



7o(ma) = -lnm-lna - / u"^" ^ ( — ^ —J— ) du 



"1 

— Inm — Ina— / v'^~^ 



. M — 1 In M 
1 1 



m(f — 1) In f 
The Corollary then follows from the factorization 

1 v^m-l k/m 



dv. (2.23) 



(2.24) 



Proposition 2. We have previously obtained the integral representation [7] (2.86) 



The idea of the proof is to suitably expand the logarithm of the integrand, and then 
to perform termwise integration. For this we write 

/ .oo .Ox eMn(l + e-^) , 

15 



- (-1)^ r £1^!:^^,+ r,-. b + Mi + e-^)] ^^ 

^ j Jo + TT^ Jo + TT^ 



= 5 + + _ g + 2S, wi - ,)|) + g ^ 

1 1 °° 1 1 °° 1 

= ^ + Ci(vr) + - E -(^ - 2Si[vr(j - 1)]) + - E "(^ " 2Si[vr(j + 1)]). (2.26) 

2 27r ^.^2 3 27r ^.^^ j 

Shifting the index in the summations gives the expression (1.20). 

Remarks. The asymptotic forms of Si and Ci are easily obtained to any order by 
repeated integration by parts. It is then easy to see that the summations in (1.20) 
have summands that are 0(l/j^), and additionally these leading terms have sign 
alternation according to (—1)-'. 

For comparison purposes, we recall an earlier result [9] 

oo 

l-7 = ^ + 2ECi(27rj), (2.27) 
^ i=i 

that readily shows how to develop 1 — 7 from 1/2 with a series of corrections, and 
the leading terms in the corrections are easily written. 
We recall that the constant (e.g., [I9], p. 345) 

1-1 = ^dt, (2.28) 

where the fractional part {t} = t—[t]. We have that Pi{x) = Bi{x—[x]) = x— [x] — 1/2, 
the first periodized Bernoulli polynomial, has the standard Fourier series pj (p. 805), 

P,(x) = -E^^^^^. (2.29) 
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Inserted into (2.22), we have 

-| oo 

= - + 25:Ci(2vrj), (2.30) 

where we integrated by parts and made a simple change of variable. 
As a second proof of (2.27) we have the following. We write 

°° foo POO foo POO p~^'"rlT 

y cosvdv e-(^+2.,>^^ _ _ / ^Qg^ / 
r"! Jo Jo Jo Jo 



-dv 



i=i 



e^'^^' - 1 

"""^^ Ul-l)- (2-31) 



/o (x2 + l)(e2'^^- 1) 2 V2 
Here the interchange of summation and integration is justfied by the absolute con- 
vergence of the x-integration. In the last step, we applied Hermite's expression for 
the digamma function ([19], p. 91 or [12], 8.361.3, p. 943). 

Further representations for combinations of the Euler constant and zeta values 
may be obtained from the following. 
Lemma 2. We have for k >2 

4.fM;.<^i-.-gMzizii. ,,32) 

We provide an operational proof, using the Dirac delta function 6. We have, 
integrating by parts. 



{ty 
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{t} 



k-1 



-dt 



\t\^-^ 1 



1 1 

3=2 3 



(2.33) 



Herein, we used d[t] = J2j ^(t ~ j)dt. 

Of course — )■ as A; — )■ oo, and we have the simple bound Ik < l/k. 

We let Bk{x) be the Bernoulli polynomials and Pk{x) = Bk{x—[x]) their periodized 
form. In applying Lemma 2 we may use the relation 



{x}" = {x- [x]y 



n + 



1 _ " /n + r 

1 r". 



Pk{x), 



(2.34) 



k=0 



with Bq{x) = 1, and the Fourier series for n > 1 [T] (p. 805) 



P .^n 2(2n-l)! - sin(2vrja:) ^ 2(2n)! ^ cos(2vrjx) 

P2n-l{X) - (-ij |,27r)2n-l fn-l ' ^Snl^j - (-ij -^^^ . 

(2.35) 



Then we have 

Corollary 9. We have (a) 

O 1 1 oo 

/2 = - - 7 - 2C(2) = ^ + E[2Ci(27rj) + jvr^ - 1 - 2j7rSi(27rj)], (2.36a) 
and (b) 

11 1 1 1 oo 

/3^^-7-5C(2)-5C(3)^j+i: 



2 1 - Trr' Ci(2ij) - - + ji^ - 2j,rSi(2,rj) 



Proof. For part (a) we use the combination 

1 ^ rcos(27rja;) sin(27rjx) 



{xY = P2{x)+P,{x)+^- = Y. 

3 = 1 



1 

3' 



(2.366) 



(2.37) 



1^ 



and carry out the integrations. For part (b) we again use the Fourier series (2.35), so 
that 



{a;}3 = P3(^)+ P2(a;) + Pi(x) + - = 



3 cos(27rja;) _j_ / 3 \^ sin(27rj2;) 



1 

(2.38) 



Proposition 3. From (e.g., [TT], p. 107) 

\nT(a) = ( a - -]\na - a + - \n(2n) - H ^-^dt, (2.39) 
V 2y 2 Jo t + a 



we have 

^(a) = lna--+/ jr^dt. (2.40) 

2a JO [t + a)'' 

Using the Fourier representation (2.29) and performing the integration gives the 
Proposition. 

Remarks. When a = 1 in (1.21) we recover (2.27), while from the case a = 1/2 
we find the companion relation 

oo 

7 + ln2 = l-2^(-l)-'Ci(7rj). (2.41) 

The sum on the right side of this equation may be determined similarly to how we 
found (2.31). 

Similar series representations in terms of Si and Ci may be obtained for Inr(a) 
and for the polygamma functions ip^^^a). For instance, from (2.39) we have 



In r(a) = ya — -j \na — a + - ln(27r) 

1 °° 1 

+— V - [2 sin(27rja)Ci(27r ja) + cos(27rja)[7r - 2Si(27rja)]] . (2.42) 
2vr . -, 7 
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From this equation at a = 1 we may conclude that 



1 OO 1 -1 

- E - 2Si(27rj)] = 1 - - ln(27r), (2.43) 
27r .^^ J 2 



and at a = 1/4 that 



,„r (j) ^ 5in(4.)-j+- E Hfl- - + 2 E + l/2)l| . 

(2.44) 

We give a second direct proof of relation (2.43). We have 

E-[vr-2Si(27rj)]=2;^-[Si(oo)-Si(27rj)]=2;^- / — 

7 = 1 J 7=1 J 7=1 J -^^""^ ^ 



-dt 



^ 1 /•°° sinf , ~ 1 /-oo , /•= 

2 y - / dt; = 2 y - / sinvdv 

fr{j Jo v + 2ti] J 7o 7o 



e 



j = i J J j = 



sint;rfw / e-^''ln(l - e-2"^)rfx = -2 / ^ ^ 

JO Jo 1 + 

°° tan~^ a; 

e^^^^ - 1 



4W„ 75^^^^- (2-45) 



In the last step we integrated by parts. We now differentiate a known result [19] (p. 
100), to write 

, 1 cos(s tan~^ t) tan^^t 

^^^^""(S-1)2+ io (l+t2)s/2 (e2-t_i)^^ 

oosin(.tan-t)ln(l + t^)^^_ (2.46) 



(l + t2)«/2 (e2-t-l) 
Since by (2.46) we have 

1 tan / 

C'(0) = --ln(2vr) = -l + 2/^ -^rft, (2.47) 
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the result (2.43) again follows. 

Proposition 4- We first demonstrate part (b), with part (a) being similar. We use 
the representation for Re x > (e.g., p3], p. 928) 

Ci(x) = 7 + In X + r i^^^l^rft. (2.48) 

JO t 



We also use the Fourier expansions 

^ cos(/3n) _ _^ 



n=l 



48 



and 



Then we have 



^ cos(/3n) 



E 

n=l 



n=l 



48 24"^ 8^^ ' 



OO -| 

e4 



n=l 



7 + ln(/3n) + 



(cost — 1) 



dt 



, 1 ^x>/.x /■Mcos(/3nt;) - 11 , 
= 7 + In /3 C 4 - C 4 + E / 

(7 + ln/3)C(4) - C'(4) + ^ (^-^t;^ + 7^(5v^ - tt^^ c/^; 
H7 + M)C(4)-C'(4)4(-^ + ^^-^). 



(2.49a) 



(2.496) 



(2.50) 



The integral over v being absolutely convergent, the interchange of summation and 
integration is justified. For (1.22d), we use the alternating zeta function for Re s > 
p. 1073, or [IS], p. 96) 



E 

n=l 
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(2.51) 



so that 



^ L_L In n = (1 - 2i-)C'(.) + 2i-(ln 2)C(s) 



(2.52) 



n=l 



For part (c) we use the summations [H] (p. 244) 



E 

n=l 



COS nx 



n 



2k 



2{2k)\ V2. 



X 



and 



E(-i) 



^cosnx (— 1) 



k-l 



n=l 



n 



2k 



2{2k)\ 



X + TT 

2n 



(2.53a) 



(2.536) 



For part (d) we use the y = case of the two summations [Hj (p. 244) 



^ cos(A;x + y) (27rY 
> ; = . csc(7ra) 



fc=l 



2r(a) 



sm[y + -)C[l-a,l-- 



na\ X 

sm I -u ] C \1 — a, — 

2 )^\ ' 27r 



(2.54a) 



and 



^, ^.fcCos(A;x + ?/) {2-kY , . 
E(-l) — = ^w7TCSc(vra) 



2T{a) 



TT ~ X 



sin(y + f)c(l-a, 2. 



7ra\ ( Ti + x 



Part (e) uses 



E 



and 



E 



; (2n + 1) 
1)" 



4-a 



(l-2-")C(a)-l, 



c(«,^)-c(a,^ 



(2n + 1)» 

together with the y = case of two summations [H] (p. 245), so that 



^ coskx (27?)" , , f /x + vra 
2^ 7^ = X csc(7ra) { — sm 



fc=0 



i2k + lY 4r(a) 



2 
22 



(2.546) 



(2.55a) 



(2.556) 



a, 



a; 

An 



+ sin 



X — TTa 



C 1 - «, 



27r — X 



X 



(2.56a) 



and 



cos/ex (27?)" f /x + vra 

— csc(7ra) <^ cos 



£^ {2k + lY 4r(a 



C 1 - a, 



47r 



C I -a, 



+ cos 



X — vra 



C I -a, 



TT — X 

47r 



C 1 - a, 



Svr — X 



Part (f) uses Pj (p. 247) 



°° 2;" cos rax _ 1 



47r 



smx 



— = 1 iz — cos x) ln(l — 2z cos x + z^) 

n[n + 1) 2z z 



■ tan 



-1 



(2.566) 



zsmx 
1 — z cos X / 

(2.57) 



Remarks. In part (a), due to the functional equation of the Riemann zeta function, 
we may write the value C'(2) = C(2)[7 + ln(27r) + 12 In A], where A ~ 1.28243 is the 
Glaisher constant, such that In A = 1/12 - C'(-l) = -[C(-l) + C'(-l)]- 

As shown by part (c), the Proposition extends to all summations I]^i(=tl)"Ci(/3ra)/'n 
with k > 1. Alternatively, successive Fourier expansions may be built up by twice 
integrating the previous one, and requiring the appropriate C{2k) constant term. 
Of course in this part, the relation (p!l], p. 14 or [12], p. 1077 or pjj p. 98) 
C(2A;) = (-l)'=-i(27r)2'=fi2fe/[2(2A;)!] applies. 

Analogous results can be given for sums J2'^=i['^ ~ 2Si{f3n)]/n'^^~^^ with k >0. 

Part (d) of the Proposition subsumes part (c), owing to the reduction C(l — "^j q) = 
—Bn{q)/n (e.g., Pj)J, p. 85) and the functional equation Bn{l — x) = (— l)"i?„(x). 
The latter relation shows the reflection property -B2fc(l — x) = B2k{x). 

Other Ci summations are possible, based for instance upon known sums for 
Er=o(±l)''«''cosfcx/(fcm + ra) [H] (pp. 246-247). 



2k 
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Summary 

We have developed a variety of series and integral representations for a family of 
functions. Prom the representations for the polygamma functions follow many special 
cases, including representations of the values ip{p/q), of harmonic numbers H^, and 
of generalized harmonic numbers H^^ . We have found representations for In F and in 
turn a product representation for the Gamma function. Most of the representations 
are quite suited for computation, and in fact their convergence may also be acceler- 
ated. The approach via the zeroth Stieltjes constant enables the development of new- 
representations, as well as the recovery of known series and integral representations 
for constants such as 7 and for the digamma function. 

We have given series representations in terms of trigonometric integrals Ci and 
Si for 7, for certain combinations of 7 and ({n)/n, for 7 + ln2, and generally for 
In r(a) and ip{a). Additionally, we have demonstrated in several ways how to perform 
summations over parameterized values of Si and Ci. 
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